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Abstract: We revise the bounds on heavy sterile neutrinos, especially in the case of their 
mixing with muon neutrinos in the charged current. We summarize the present experimen- 
tal limits, and we reanalyze the existing data from the accelerator neutrino experiments 
and from Super-Kamiokande to set new bounds on a heavy sterile neutrino in the range 
of masses from 8 MeV to 390 MeV. We also discuss how the future accelerator neutrino 
experiments can improve the present limits. 
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1. Introduction 

The existence of three active neutrinos and of the mass eigenstates v\ , z^, ^3 is well es- 
tablished, but the existence of a forth heavy mass eigenstate, v^, mainly in the direction 
of a sterile neutrino, remains an open question. Sterile neutrinos are SU(3)xSU(2)xU(l) 
singlet fermions, which can mix with ordinary neutrinos. 

There many reasons why one is interested in the limits on sterile neutrinos. Heavy 
mostly-sterile neutrinos have been investigated for their role in cosmology and astrophysics 
(see, e.g., Ref. Q). A keV sterile neutrino is a viable dark matter candidate [Q, which can 
also explain the origin of the pulsar kicks 0] . Sterile neutrinos provide a viable framework 
for baryogenesis ||, |g]. It was pointed out |zj that if such heavy neutrinos constitute a small 
but non-negligible fraction of dark matter, their decays into e + -e~ pairs might produce the 
511 keV gamma line observed by the INTEGRAL 7-ray observatory. Decays of heavier 
neutrinos have been proposed to explain the early ionization of the Universe [^ . If neutrinos 
are Majorana particles, they could mediate processes that violate the lepton number by 
two units, such as neutrinoless double beta decay, muon-positron conversion, and rare kaon 
decays §. In particular, for masses 245 MeV < < 388 MeV, the decay K + —* 7r~/i + / u + 
could be strongly enhanced pQ]. From the theoretical point of view, heavy sterile neutrinos 



with masses in this range arise naturally in extended technicolor models [11]. Because of 
these interesting possibilities, we want to map out the parameter space available for sterile 
neutrino masses and mixing parameters. 

Here we consider the bounds on heavy neutrinos with masses 1 MeV ^ ^ 400 MeV, 
which are produced in pion, muon and kaon decays. We will review the existing limits on 



the mass and mixing of a heavy, mostly-sterile, neutrino [12, 12, 14, 15, 16, 17, 18, IS, 20 
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21, 22, 23, 24, 25, 26, 27 1. We will then derive new bounds based on a re-analyses of the 



data from neutrino oscillations and accelerator neutrino experiments. In the 10-390 MeV 



mass range we use the data from accelerator experiments p4] , p5| , p6| , 27 1. For masses 
8 MeV < m h < 105 MeV, we use the Super-Kamiokande (Super-K) data |2§]. We will 
also discuss ways in which the present and future neutrino oscillation experiments, such as 
MiniBOONE pi, K2K p|] and MINOS plj, can strengthen the present bounds. We also 



review the best current limits, based on big-bang nucleosynthesis (BBN) |32|j . 

Let us characterize the mixing of the heavy neutrinos with the active neutrinos v a 
(a = e, fi, t) by the corresponding element in the mixing matrix U, which is the mixing ma- 
trix between the electroweak eigenstates and the mass eigenstates. The matrix U accounts 
for mixing in the neutral current (NC) interactions. In processes mediated by the lepton 
charged current (CC), the matrix U enters in combination with a unitary matrix V, which 
diagonalizes the charged lepton mass matrix. In principle, CC and NC interactions allow 
one to measure or constrain separately the elements of VU and U, respectively. If one 
takes V = 1, the elements of the matrix U can be interpreted as neutrino mixing angles. 
We will not make any assumptions about the matrix V and will present the constraints 
in full generality. We will assume, for simplicity, throughout our analysis that mixes 
mainly with v^, in the charged current, (VU) e h, (VU) T h ~ 0, while we allow for U a h ^ 0, 
a = e,fJL, t, in the neutral current. 

In Section 2 we review the bounds coming from the analysis of the spectrum of muons in 
pion and kaon decays. In Section 3 we study the case of mostly-sterile neutrino decays into 
visible decay products which would be observable in a detector. We discuss the limits on 
the mixing which can be obtained from accelerator neutrino and Superkamiokande data and 
the prospects for strengthening such bounds in present and future neutrino experiments. 
In Section 4 we briefly review the bounds from big bang nucleosynthesis and supernovae. 
Finally, we summarize our results in the Conclusions. 

2. Bounds on production 

The different massive neutrinos z/j, i = 1,2,3,4, if they exist, are produced in meson 
decays, e.g. tt^ — > (J, i>i, with probabilities that depend on the mixing in the charged 
current, VU. In our analysis we will assume that the heavy sterile neutrinos mix mainly 
with muon neutrinos in the charged current. We allow for mixing with all active neutrinos 
in the neutral current. The energy spectrum of muons in such decays would contain 



monochromatic lines Jl2j at 

Ti = (ml + - 2m^m^ - ml i )/2m 7r , (2.1) 

as long as Tj > 0. Here Tj is the muon kinetic energy; m n , m^, and m„. are the masses 
of pion, muon and the ith neutrino mass eigenstate, respectively. The dominant line is 
obtained for nearly massless neutrinos, ^1,2,3; at Tq = 4.120 MeV. Additional peaks in 
the muon energy spectrum would be present at a position related to the mass of the heavy 
neutrino, Eq. ( |2.lD , and with a branching ratio that depends on the mixing angle. The 



same is true for muons from K decays. Searches for peaks in pion decays [12, 13, 14, 15 
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Figure 1: The exclusion plot for |(W) M /j| 2 based on the energy spectrum of muons in pion decays. 
The excluded region is indicated in gray color (magenta color online). The bounds are taken from 
the analysis reported i) in Ref.JIJ] for the dash-dotted line indicated as "[1]"; ii) in Ref.Jll| for the 
dotted line "[2]"; iii) in Ref. |l6| for the solid line shown as "[3]"; iv) in Ref. Jl7| for the dashed line 
"[4]"; v) in Ref.Q for the dashcd-double dotted line labelled "[5]". The bounds are 90% C.L., 
except for the one marked "[5]", which is 95% C.L. 



0, 0, |T|, H, H and in kaon decays ||, || found no signal and set stringent bounds 
on \iyU) iJL } l \ 2 for masses ^ 360 MeV. The corresponding excluded regions from pion 
and kaon decays are shown in Figs. 1 and 2 as gray regions, respectively. Different lines 
represent limits obtained in different experiments (see captions of Figs. 1 and 2). 



3. Bounds on decays 

A heavy neutrino produced in tt^ and decays would subsequently decay, and its decay 
products could be detected. The absence of such detection translates into strong bounds 
on the mixing angles with active neutrinos. In the mass range 10-390 MeV, we set a new 
bound on \{yU)^hU a h\-, as shown in Fig. ||. 

A heavy neutrino mixed with active neutrinos can decay into different channels de- 
pending on its mass. If 1 MeV ^ m/j ^ 105 MeV, Vh decays via neutral currents into (i) 
an active neutrino and an electron-positron pair, Vh — > v a + e + + e~ , (visible channel) and 
(ii) into three neutrinos, — ► u a + u.j + Pj, (invisible channel). 

In the simplest case, when sterile neutrinos directly couple only to one active neutrino 
in the neutral current with mixing U a h , the mass of the heavy neutrino is much larger 

1 In the literature the mixing is often parametrized by a mixing angle 8. We have that sin 2 8 = |[/ a h| 2 - 
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Figure 2: The same as in Fig. [j] but for i>h from kaon decays. The bounds are taken from Ref. |2^] 
for the dash dotted line indicated as " [1]" at 2 a and from the data of the experiment by R. S. Hayano 
et al, Ref. §|, for the solid linc labelled "[2]" at 90% C.L. 



then electron mass m e <C m/j < m^, and we assume that in the CC current (VU) e h ~ 0, 
the decay width is given by [33]: 



r _ 1 + (g£) 2 + 9R n 2 _5 \jr |2 /o-i\ 
h ~ 7687T3 (^F m h\ U ah\ , (3.1) 

where g e L = 1/2 + sin 2 ^, g^ T = — 1/2 + sin 2 ^, g# = sm 2 9-^. If the heavy neutrino 
mixes with v e in the charged current, the decay — > v e e + e~ would receive an additional 
contribution from CC interations. In the absence of charged currents, the corresponding 
branching ratios are equal to B e = {g — l)/g ~ 0.11 for the visible channel and B v = 
1/g ~ 0.89, for the invisible channel. Here we have defined g = 1 + (gl) 2 + g\- For 
larger masses, new channels are open, namely Vh - > and Uh — * n + l±~v a which are 

mediated by neutral currents, and — > fj,~e + u e , — > fi + pTv^ and — ► 7r + /i~, due to 
CC interactions. For > 77%, the two-body decays dominate and typically the halfdife 
time is of order r h ~ (10~ 8 - 10~ 9 ) |t/^|- 2 (|(V£/) M/l |- 2 )s, for the NC (CC) mediated 
processes. 

3.1 Searches for accelerator neutrino decays 

The decay length of a sterile neutrino with energy S> is = c r^jp, where 
"fp = Eh/mh is the Vh gamma factor. The fraction of heavy neutrinos that can reach the 
detector before decaying is exp (— R cr /L^), where R CT is the distance from the neutrino 
production site to the detector. Of these neutrinos, for h/L^ <C 1, a fraction B^h/L^ 
decays in the detector via a visible channel. Here B v [ s is the branching ratio of the given 
decay channel and h is the length of the detector. The number of heavy neutrino decays 
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Figure 3: The bounds on \{VU) fi h(VU) e h\ versus mh obtained from searches of Vh decays. The 
excluded region is indicated in gray color (magenta color online). The limits are taken for the lines 
labelled as i) " [1]" (dashed-dotted line) from the data in Ref. |J] ; ii) " [2]" (dashed-double dotted 
line) from the experiment in Ref. (27); iii) "[3]" (solid line) from Refs. @ |§; iv) "[4]" and "[5]" 
(dashed lines) from the analysis in Ref. p(|; v) "[6]" (dotted line) from the experiment reported in 
Ref. pi. The limits are at 90% C.L.. 



in the detector is then given by (see, e.g., Ref. [ 

N = N 1V ^ K B(M 



/j,v h )B vis —tte 



(3.2) 



where k is the number of pions and kaons, B(M — * \iv^) is the branching ratio of the 
meson decays into a muon and a heavy neutrino, f2 and e are the detector acceptance and 
efficiency, respectively. In a similar way one can consider heavy neutrinos produced in 
muon decays. 

Limits on the mixing of the heavy neutrino with v e and v„ were set by different 
]. We review them in Fig. y. In the mass range 1 MeV < < 



experiments |24|, [25], |26|, 
33.9 MeV the excluded region comes from heavy neutrino production in pion decays. For 
higher masses, 40 MeV ^ ^ 360 MeV, kaon decays were taken into account. Different 
visible decay channels have been studied. The channel K + — > fi + Vh M + (l JL ~ e+ u e) + c.c. 
was used to constrain the elements of the mixing matrix VU, for masses up to ~260 MeV 
(see Fig. ||). For heavier masses a new decay channel is open, K + — > n + Vh — > ^ + (7r + /i _ ) + 
c.c, and dominates. 



We note that in Refs. p4| , 25, 26, 27] the NC contribution to the decay of Vh has not 
been taken into account. They mediate the principal decay modes both for neutrinos with 
rrih < m n o: (y^ — > vvv) and with nih > m^o: [y^ — > it® v.) However this omission does not 
affect the bounds obtained in Refs. [24, |25|, |26|, ||J 



because their visible decay channel is 
always dominated by the charged current (CC) interactions. Since the flavor of the final 
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Figure 4: The same as in Fig. || but for the bounds on \{VU) ll ,hU a h\- For the lines labelled as " [1]"- 
"[4]", the limits are obtained from a reanalysis of the data reported in the references as indicated 
for Fig. |. 



neutrino is not detected in these experiments, the channels K + —* — * /i + (e - e + v e ) , 
mediated by CC, and K + — > p^v^ — > p + {e~e + v e ^^ T ), mediated by NC, are effectively 
indistinguishable. The former channel was used to constrain the \(VU) e h(VU) ^ mixing. 
However, the same data can be used for setting a bound on \(VU)fj,hU a h\, with a = e,/j,,T, 
if one includes the contribution due to the latter channels. 



We used the limit on \{VU) eh {VU)^ lh \ from Refs. |2|, ||, g |7||, corrected it by 
including the contribution from the NC to the total decay width, with (VU) e h negligible, 
and the branching ratio of K + — > n Vh ~^ P + ( e ~s + i' a ) channel, and translated it into 
the new bound on \(VU) fl f l U a h\- These bounds have not been previously discussed in the 
literature and are reported in Fig. |I[ The new limits are given at the same confidence level 
as the ones on \(VU) e h(VU)^h\. 

If V = 1, our limit on \U^h\ 2 turns out to be the strongest limit in the mass range 
34 MeV < m h < 200 MeV (see Fig. |). 

Analogously, the bounds on | CTe/i. | 2 [25, 26| from ir + ,K + — > e + Vh — > e + {e + e~v e ) + 
c.c, mediated by CC interactions, can be used to constrain \(VU) e hU a h\ from the decays 
ir + ,K + — > e + Vh — > e + {e + e~v a ) + c.c, which are induced by neutral currents. Also these 
bounds are new but a detailed analysis is beyond the scope of the present article. 

3.2 Decays of atmospheric heavy neutrinos 

We point out that one can also set an independent limit based on non-observation of 
atmospheric sterile neutrino decays by Super-Kamiokande. In the following we provide a 
qualitative analysis of Super-Kamiokande data in order to put a bound on the decays of 
heavy neutrinos, if produced in the atmosphere. We also give an estimate for the limit on 
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Figure 5: The strongest bound (indicated as line "[3]" in Fig. 4 from the reanalysis of data from 
Ref. [^5[ |2(| (diagonally hatched region with dashed-dotted contours). The limits from big bang 
nucleosynthesis are indicated with the light blue-light gray region with continuous contours, and the 
previous bounds from experimental searches are shown as dark gray regions with dashed contours. 
For simplicity, here we take VU = U. 



\U a h\ 2 , which we expect to be correct within a factor of 2-3. We assume for simplicity that 
(VU)^h dominates in the charged current while (VU) e h ~ 0. 

Heavy neutrinos can be copiously produced in the atmosphere and can decay inside the 
Super-K detector generating such a high event rate. The differential flux dF^/dE of heavy 
sterile neutrinos produced in the atmosphere in pion, kaon and muon decays is related to 
the active neutrino flux F a (a = e,fi,r): 

Mm. = (3. 3 ) 

where dF a (E) / dE is the differential flux of (active) atmospheric neutrinos. 
The expected rate of decays, R, detected in a given energy bin is 

R(J*0 = B e e{E x ) | j E d E^A^-e-^ j , (3.4) 

where B e is the branching ratio of the — * e + e~v a channel, due to NC, and e(-Ei) is the 
efficiency to detect the decay products visible in the detector, E\ is their average energy, 
A is the detector area. This rate should not exceed the rate of events observed by the 
Super-K experiment per energy bin, R(-Ei) < -Revent C^i ) ■ 

It is clear that the highest number of events should come from neutrinos with the 
smaller gamma-factor , ~ 2 — 3 (we conservatively take jf ~ 2 to ensure that the 
production of is not suppressed by threshold effects). If the gamma factor of is small, 
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Figure 6: The exclusion plot based on the non-observation of e-like events from the heavy Dirac 
neutrino decays in Super-Kamiokande detector at a rate higher than the observed 7r°-like event 
rate (region with red dash-dotted contour). Also shown are the bounds from big bang nucle- 
osynthesis (light blue- light gray region with continuous contours), and from experimental searches 
(regions with dashed contours). We assume VU = U. 



"fF ~ 2 — 4, most of the emitted e~ and e + produce in the Super-K detector 2 separate 
Cerenkov-light cones, which would be interpreted as 2 e-like events coming from the decay 
of 7T° — > 77 (7T°-like events). We use the data reported in Ref. [28] where the invariant 
mass distribution for 7r°-like events is shown per energy bin of 10 MeV, for 1489.2 days 
and 22.5 kton of fiducial volume, V ~ Ah. In the case of a three-body decay, the invariant 
mass reconstructed from the energy and momentum of e + and e~ is 



m\ 2 



(Pi + Vn) 



2 i 2 

m h + m u 



2m h E v 



(3.5) 



where are the four-momenta of e + and e , and m v <C and E v are the mass and 
the energy carried away by the undetected in the reference frame of v^. We assume, 



for simplicity, that the heavy neutrino decays in 3 relativistic particles, 



e + and e 



each having, on average, 1/3 of the total energy. The reconstructed invariant mass can be 
related to the heavy neutrino mass m\ 2 — 1/3 m^. 

We compute the number of 2e-like events in each bin of mf 2 expected from Vh decays, 
with ~ 2 — 4. We take into account the energy dependence of dFh/dE: the differential 
flux is with good approximation constant for E ^ 100 MeV, and then it decreases (see, e.g., 
Ref. J34}|). Comparing the number of expected events from the decays with the 7r°-like 



events, -/V exp , in each invariant mass m\ 2 bi n > we obtain an upper bound on | {VU) 



ah I • 



We show this bound in Fig. [6| (vertically hatched region), if V = 1, for Dirac sterile 
neutrinos which mix mainly with Vn- For Majorana neutrinos the bound is somewhat 
stronger because the decay rate is larger, as they could decay both in one channel and in 
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its CP-conjugate. A more detailed analysis which includes the contribution of decays of 
with < 2, threshold effects in the production of Vh in pion and muon decays, the zenith 
angle dependence and 7i?-dependence of the suppression factor exp(—R cr / La), should be 
performed but it is not the main focus of our paper. However we expect our conservative 
bound to remain valid within a factor of a few, as it would be strengthened by the inclusion 
of events with 7^ < 2, the suppression exp(— R cr /Ld) is not very important for the masses 
and mixing parameters considered, the number of decays with large 7_f which could be 
misinterpreted as single e-like events is suppressed by a larger and, for large masses, by 
a smaller flux of Uh- 

Most of the decay products e + e~ from heavy neutrinos with larger gamma factors, 
7f 3> 3, would produce 2 nearly overlapping Cerenkov rings, which would be recorded 
as a single edike event with twice the energy. Limits on the allowed range of rrih and 
\(VU)fj,hU a h\ based on these events are typically of the same order of magnitude but are 
somewhat weaker than those shown in Fig. |6|. 

For masses ^7 — 8 MeV, the Super-K threshold limits one's ability to set a bound 
because only neutrinos with Eh ^> rrih would be above the threshold. 

In principle one could use the same technique to constrain the mixing angle for heavy 
sterile neutrino with masses up to mh < tuk — m^, produced in K decays. However, this 
bound would be very weak because too few neutrinos come from K decays. Our bound 
based on the Super-K data is a factor of few weaker than the bound that we obtained 
by re-analyzing the accelerator data, as discussed above. A much tighter bound, possibly 
stronger than the present limit for masses in the 8-105 MeV range, could probably be 
obtained if one included the directional and other information from the Super-K data not 
available to us. The bound can be further improved in the future, after more data is 
collected by Super-K or by a bigger detector, e. g., Hyper-K. 

3.3 Future accelerator experiments 

Present and future accelerator neutrino experiments KEK to Kamioka (K2K) p0|, Mini- 



B00NE [^9| and MINOS |31], have the possibility of setting new, possibly stronger, bounds 
on the mixing parameters. Neutrinos are produced in the decay of pions, muons and kaons. 
The number of proton on target is very high: 10 20 for K2K and 10 21 for MiniBooNE and 
MINOS. The average neutrino energies are around 1 GeV for K2K and MiniBooNE and 
2.5-3 GeV for MINOS. If heavy neutrinos are produced, they would travel toward the 
detectors where the e^, // ± , 7r° and 7r ± generated in their decays would be detected. In 



MiniBooNE 2£], the detector is of length ~ 10 m, and collects the Cerenkov and scintilla- 
tion lights. Both K2K and MINOS have a near and a far detector which can both be used 
to search for heavy neutrino decays. The K2K long-baseline experiment [ 30 1 uses a near 



1 Kton water Cherenkov detector and as a far detector Super-K. The MINOS experi- 
ment has a near and a far iron-scintillator detectors of length h ~ 20 — 30 m, respectively. 
For the far detector timing can be used to reduce the backgrounds. The neutrino beam is 
pulsed with a typical duration of the spill of few /is. The data read-out of the signal in the 
detector is done in a time-window around the beam spill of 1.5 and 10-20 fis for K2K and 
MINOS, respectively. The heavy neutrinos will take longer to reach the detector, typically 



-9- 



few to tens of ^s after the light neutrinos. Therefore a sizable fraction, if not most, of the 
heavy neutrino decay signal could be time-separated from the background given by the 
accelerator v e and interactions. Directionality can be used to disentangle a signature 
from active atmospheric neutrinos. 

As discussed previously, depending on their masses and on the mixing with the active 
neutrinos, the heavy neutrinos decay mainly in — > 3f a , fh — ► e + e~v a , t>h — > ^e^u^, 
and for rrih > m w in 7r°z/ a , ir^e^, it^/j^. The CC induced decays would allow to constrain 
the product \(VU) e h(VU)nh\, while the NC processes \{VU') tl hU a h\- The flux contains 
with different 7_f, which would give different signatures in the detector. A detailed analysis 
which takes into account the heavy neutrino spectrum, the different decay channels and 
their signals in the detector, the relative backgrounds needs to be performed in order to 
achieve a good evaluation of the sensitivities for these experiments. Nevertheless we can 
provide an order of magnitude estimate for the limits which can be reached on the mixing 
parameters. Searching for 2 e-like, single e-like, /x-like and pion events, K2K, MiniBooNE 
and MINOS should be able to reach sensitivities as good as a fewx 10~ 7 for rrih ~ 100 MeV. 
The sensitivity would be of the order of a fewxlO -9 and a fewxl0~ 10 for heavy neutrino 
masses rrih ~ 200 MeV and 300 MeV, respectively. 

In contrast with the other experiments, MINOS has a very good discrimination of /i + 
and Neutrinos with rrih ~ 210 MeV would decay in this channel with a branching 
ratio which is typically ~ 10~ 3 . The decay Vh — > /i + ' yTv^ receives a contribution from CC 
interactions and therefore allow to constrain not only \(VU)^hU a h\ but also \(VU)^h\ 2 - As 
no significant background is expected for this decay mode, this would be a clear signature 
of the existence of heavy neutrinos. The number of events would allow to establish the 
value of the relevant mixing angle and the reconstructed invariant mass the value of m^. 

We note that a possible future detection of the peak in the invariant mass distribution 
for 2 e-like events from — ► e + e~v a , after the contribution due to 7r°-decays is sub- 
tracted, for 2 /[/-like events due to Vh P" V a and for 7r ± and e T produced in 
the decays Vh — * 7r ± e T ( / u =F ) decays, would signal the existence of a sterile neutrino and 
would allow a measurement of its mass and mixing. If the different decay channels can 
be well distinguished, a comparison of the observed branching ratios with the theoreti- 
cally predicted ones would allow to confirm the hypothesis of heavy neutrino decays as 
the origin of the observed signal. Furthermore, if from indipendent measurements we have 
that \{y"U) ll hU tl h\ 7^ liVU)^] 2 , it would be possible to establish that V ^ 1 and test the 
existence of mixing in the charged lepton sector. 

As already noticed in the previous discussion, the data from the MiniBooNE and MI- 
NOS experiments can be used to constrain also the combinations |(l A J7) e / l | 2 , |(V?7) e ft(VC/) /J / l 
\(VU) e hU a h\, once the proper decay channels and branching ratios are taken into account. 
Also for these mixing terms, we expect a strengthening of the present bounds. 

4. Bounds from cosmology and astrophysics 

Heavy sterile neutrinos in the MeV mass range would be produced in the Early Universe 
and subsequently decay affecting the predictions of Big-Bang Nucleosynthesis (BBN) for 
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the abundance of light elements and in particular of 4 He (see, e.g., Ref. [p2[ ). The main 
effect would be to increase the energy density, leading to a faster expansion of the Universe 
and to a ealier freeze out of the n/p-ratio. In addition, the decay of Vh into light neutrinos, 
in particular, v e , would modify their spectrum and the equilibrium of the n — p reactions. 
We report in Figs. 5 and 6, the bounds on the masses and mixing angle of the heavy sterile 
neutrino, which can be obtained from big bang nucleosynthesis. Note that in Ref. [32] the 
BBN bounds were derived up to nih 200 MeV. It was shown that in the region 140-200 
MeV the dominant decay chanel is into pions. We extrapolate those bounds above the 
QCD phase transition up to 400 MeV (dashed border of the BBN excluded region) taking 
into account the change in relativistic degrees of freedom, which softens the BBN bounds. 
In order to produce a more careful BBN bound for > 200 MeV a detailed analysis 
similar to the one in Ref. [32] should be performed but is beyond the scope of the present 
study. In Figs. || and|6] we indicate our estimated BBN bound above 200 MeV with dashed 
lines to underline that these limits should be considered valid within factors of a few. 

In principle, SN1987A could be used to exclude sterile neutrinos with mixing angles 
10 - ' ~ \U^h\ 2 ~ 10~ 2 and masses m s ~ T com , where T core = 30 — 80 MeV is the core 
temperature of the neutron star at about 0.1 second after the onset of the supernova 
explosion [32|. For masses in excess of T core , the production of sterile neutrinos is suppressed 
by the Boltzmann factor. The emission of sterile neutrinos from the core depends on the 
pattern of their mixing with active neutrinos, and in some cases the emission history can be 
very complicated [pq ]. Given the uncertainty in the SN models, one cannot set a reliable 
bound outside the range already excluded by the combined BBN and new accelerator 
neutrino bounds shown in Fig. ||. 



5. Conclusions 



We have considered heavy sterile neutrinos mixed mainly with muon neutrinos in the 
charged current. These neutrinos would be produced in pion and kaon decays due to 
their mixing with f„. They would subsequently decay into standard model particles, e.g. 
neutrinos, electrons and positrons. The non-observation of the decay products in past 
dedicated experiments has sets some stringent bounds on the mixing angle between heavy 
and active neutrinos. We have reanalyzed the existing accelerator data and set new bounds 
on heavy neutrinos with masses in the range 10 MeV ^ mn ~ 390 MeV. In addition, we 
have used the Super-Kamiokande data to set a new independent bound on the mixing of a 
sterile neutrino with mass in the 8-105 MeV range. We have also discussed the potential 
of future experiments, as K2K, MiniBooNE and MINOS, in improving the present limits. 
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